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ABSTRACT 

We investigate the motion of the primary beam outside the light cylinder in the pulsar 
wind. Inside the light cylinder both primary and secondary plasma move along dipole 
magnetic field lines where their energies can be arbitrary. But at larger distances the 
theory predicts quasi-radial motion with the velocity exactly corresponding to the 
drift velocity which cannot be the same for primary and secondary plasma. Hence, 
the deceleration of the primary beam is to take place simultaneously resulting in the 
acceleration of the secondary plasma. We investigate this process in the three-fluid 
MHD approximation and demonstrate that for most pulsars the energy of the beam 
remains practically unchanged. Only for young radio pulsars (Crab, Vela) essential 
deceleration up to the energy of the secondary plasma takes place outside the fast 
magnetosonic surface rp ~ (10-100)i?L, the energy of secondary plasma itself increas¬ 
ing insufficiently. 
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1 INTRODUCTION 

According to the common point of view, the activity of ra¬ 
dio pulsars is connected with the electron-positron plasma 
generated near magnetic poles (Manchester & Taylor 1977; 
Smith 1977). Numerous works devoted to this subject (Stur- 
rock 1971; Ruderman & Sutherland 1975; Arons 1981; 
Daugherty & Harding 1982; Gurevich & Istomin 1985; Is¬ 
tomin & Sobyanin 2007; Medin & Lai 2007; Timokhin & 
Arons 2013) demonstrate that the outflowing plasma con¬ 
sists of the primary beam with the Lorentz-factor 7 b ~ 10 ' 
and the number density n b « ugj, where ugj = £lB/2irce 
is the Goldreich-Julian number density, and the secondary 
electron-positron plasma with nr ~ (10 3 -10 5 )nGj, and 
7 ± - 10 2 . 

It is clear that inside the light cylinder r < Rl = c/Q 
both the primary beam and the secondary plasma move 
along the dipole magnetic field lines where their energies 
can be arbitrary. On the other hand, at large enough dis¬ 
tances from a neutron star r Rl both the analyti¬ 
cal theory (Michel 1994; Beskin et al. 1998) and numeri¬ 
cal simulations (Contopoulos et al. 1999; Spitkovsky 2006; 
Tchekhovskoy et al. 2013; Philippov & Spitkovsky 2014) 
predict the quasi-radial outflow of the relativistic electron- 
positron plasma along the poloidal magnetic field with the 
radial velocity exactly corresponding to the drift motion 
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Clearly, this drift velocity cannot describe simultaneously 
the motion to the beam and the secondary plasma which 
have sufficiently different energies. 

Thus, the question arises about the energy evolution 
of these two components as they escape the pulsar magne¬ 
tosphere and propagate outwards forming the pulsar wind. 
We show that efficient deceleration of primary beam can 
be realized for young energetic radio pulsars (Crab, Vela) 
only. In this case the deceleration up to the energy of the 
secondary plasma takes place outside the fast magnetosonic 
surface tf ~ (10-100) Rl, the energy of secondary plasma it¬ 
self increases insufficiently. For ordinary pulsars with period 
R ~ Is, the energy of the primary beam (due to escaping 
to the region with rather small electric and magnetic fields) 
remains practically unchanged. 

The paper is organized as follows. We start with the 
description of basic equations of three-fluid MHD in Sect. 2. 
In Sect. 3 we describe the general properties of these equa¬ 
tions. The beam damping is discussed in Sect. 4. Finally, in 
Sect. 5 the main results of our consideration are formulated. 


2 BASIC EQUATIONS 

In this section we derive the basic three-fluid MHD equa¬ 
tions. Throughout the paper we use the spherical coordi¬ 
nate system (r, 9, ip) with unit vectors (e r , eg, e v ). All quan¬ 
tities having the superscript ’b’ correspond to the primary 
beam, and the quantities with superscript ’±’ refer to the 
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secondary plasma particles (electrons or positrons). Finally, 
electron charge equals —e. 

To describe the deceleration of the primary beam (and 
possible acceleration of the secondary plasma) we start from 
the well-known Michel (1994) monopole force-free solution 


B r = B h ^, 

Eg = B v = — Bl— ^ sin#, 
r 


( 2 ) 

( 3 ) 


corresponding to the zero particle mass m e . Within this ap¬ 
proximation, the particles move radially with the velocity 
v = c, so that the current density j r = cp e and the charge 
density 


Pe 


Q,B r 

27rc 


cos # 


( 4 ) 


just corresponding to Goldreich-Julian charge density 


S7B 


Pgj = — 


27TC 


( 5 ) 


One can easily check that the fields (2)-(3) with the 
charge density p e (4) and the current density j r = cp e are the 
exact solutions of the time-independent Maxwell equations 


V ■ E = 47rp e , V x E = 0, 

(6) 

V • B = 0, V x B = — j. 

c 

(7) 

For this reason one can model the outflow with the massless 

particles moving radially with the velocity v 
having the following number density 

= ce r and 

„± i Ml Rl 

2nce r 2 ’ 

(8) 

b H-Bl Rl n 

n = -- cos 0. 

Z7rce r z 

(9) 


Here A = en e /|pGj| i§> 1 is the so-called multiplication pa¬ 
rameter which equals 10 3 -10 5 for most radio pulsars. In 
what follows we consider for simplicity the case of A = const. 
Such a choice corresponds to the equality of the secondary 
electron and positron number densities. 

As one can see, the only difference with our previous 
works (Beskin & Rafikov 2000; Beskin et al. 2004) is that we 
consider here not two-component (i.e., secondary electron- 
positron plasma) but three-component outflow containing 
also the primary beam having Goldreich-Julian number den¬ 
sity n h = | pgj |/e and Lorentz-factor 7 b ~ 10 7 3> 7 *. As is 
shown in Appendix A, this value is in agreement with the ra¬ 
diation reaction force acting on the beam particle as it moves 
along dipole magnetic field line inside the light cylinder. On 
the other hand, as is shown in Appendix B, two-stream in¬ 
stability is not effective, and, hence, three-fluid hydrody- 
namical approximation under consideration is good enough 
to describe the evolution of the primary beam. 

Then, following the papers mentioned above we can in¬ 
clude into consideration the finite particle mass as a small 
disturbance of the force-free solution. One can do it because 
for magnetically dominated outflow a 1, the outflow re¬ 
mains actually radial (Tomimatsu 1994; Beskin et al. 1998). 
Here 


fleB^Rl 

4A me 3 


( 10 ) 


is the Michel magnetization parameter. For ordinary pulsars 


(P ~ 1 s, .Bo ~ 10 12 G), one has a ~ 10 3 — 10 4 , and only for 
the fast ones (P ~ 0.1 — 0.01 s, Bo ~ 10 13 G), a ~ 10 5 — 10 6 . 

As a result, to find the structure of the flow we have 
to analyze Maxwell equations (6)—(7) and the equation of 
motion for all three components: 

(v^VJp* = ±e(^E + ^ xB^) , (11) 

(v b , V)p b = -e x B^ . (12) 


It is convenient to introduce electric potential < f? e (r, 6) and 
the magnetic flux function X&(r, #), so that 


E= — V4>eM), 

B ^ = 

p 2ivr sin # 

Then, one can seek the solution in the form 


(13) 

(14) 


nA 

VLBl i?L r \ i ± ( 

= 2ncer 2 ^ + V ^’ 

(15) 

n h 

HBl Rl r „ , b , m l 

= 2ncer 2 [ mSe + V 

(16) 

Mr, 9) 

HRlBl r a , xr mi 

= -[— cos 9 + o(r, 0)\ , 

(17) 

9) 

= 27tRlRl [1 — cos# + e/(r,#)] , 

(18) 

±,b 

v r 

= c[l-^’ b (r,#)] , 

(19) 

±,b 

v o,<p 

= <^ b (r,#), 

(20) 


with the small disturbances r/, 5, ef and £ corresponding to 
the small perturbations caused by the finite particle mass. 
Accordingly, using (13) and (14), one can write down 


B r 

Bg 

Bcp 

E r 

Eg 


Rl 


= Bl 


= -eB L 


+ 

Rl 


df N 


rsin# dr 
Rl 

3L- 


sin 9 89 / 

df 


[-sin# - C(»",#)] , 


o Ml OS 

~-DL- TT-, 


( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 


Here all deflecting functions are supposed to be <§; 1. For 
v r = c, vg = 0, and v v = 0 for all types of particle, we 
return to the force-free Michel solution. 

Now, substituting Eqns. (15)-(25) into Maxwell 
equtions (6)-(7), we obtain to the first-order approximation 
the following system of linear equations: 


b N d ( 2 d5\ 
2 (A7 + r ? ) + -^ -j + 

d^(Csin#) = AAe.- 

= \ (AACe-Cecos #) , 

e 8 2 f e 8 ( 1 df 
sin # dr 2 r 2 89 \ sin # 89 


1 d ( ■ a 05 ' 
sin 9 86 \ Sm ° 89 / 

£ b cos # — Ap — rj h 


= 0, 




(26) 

(27) 

(28) 

(29) 


Accordingly, equations of motion (11) and (12) now looks 
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like 


d~r ^ 7 ) + ^- = 

= ± 4A.(-I§ + <- 


sin# ± c ± 
r ?r + flr 2 ^ 


(30) 


_a 

dr 


(* b ) 


«.V 


I 1 ^ , C s in # | c ^ b N 

a rd9 + r r ^ r + flr 2 ^ 3 1 ’ 


d_ 

dr 

d 


sM 


w / ±x ... ( dS sin 6 ± 

- ( 7 ‘) = ±4A„ ( -- - — 


= —4Act l -jr- — ———£<? 

or r 




= =f4A<t e 


df 


fir sin# dr 


" fir 2 ^ ) ’ 


& (bA 
d~r V - 7 ) 


+ 


&l h 


4AfJ * £ fir sin# dr + flr 2 ^ ^ ' 


(31) 

(32) 

(33) 

(34) 

(35) 


Here A A = A + — A ~, so these terms vanish in one fluid 
MHD description. 

Formally, the system of equations (26)-(35) requires fif¬ 
teen boundary conditions. We consider for simplicity the 
case f IR/c -C 1 when the star radius R is much smaller than 
the light cylinder. As a result, one can write down the first 
nine boundary conditions as 


£(Rl,0) = 0 , (36) 

&(RuO) = 0, (37) 

O-Rl,#) = O’ (38) 

£(Al,#) = 0, (39) 

7 ± (i?L,#) = Tin, (40) 

7 b (#L,#) = Tin, (41) 


where Tin and J;/ are the initial Lorentz-factors of secondary 
plasma and beam particles correspondingly. 

Finally, as a boundary conditions one can put 


<5(Rl, #) 

= 0 , 

(42) 

ef(Rh,0) 

= 0, 

(43) 

r]- (Rl, 8) 

= 0, 

(44) 

V b (^L,#) 

= 0, 

(45) 


resulting from the relation El + (fli?L/c)e v x Bl = 0 which 
implies a rigid rotation and a perfect conductivity of the 
surface of a star. It is necessary to stress that all the ex¬ 
pressions are valid only for a quasi-monopole outflow which 
exists only outside the light cylinder. 


3 GENERAL PROPERTIES 

It is clear that in general case the analytical solution of Eqns. 
(26)-(35) cannot be found. On the other hand, this system 
of equation contains some integrals of motion which allow us 
to obtain the necessary information about the energy of the 


outflowing particles. In particular, as was demonstrate by 
Beskin & Rafikov (2000), in the two-fluid approximation this 
approach reproduces well-known one-fluid MHD results for 
the position of the fast magnetosonic surface for T; n -C <r 1,/3 

tf = ct 1 ^ 3 sin ” 1,/3 # Rl (46) 

and the Lorentz-factor at this surface tf = T^J’f) 


Tf = a 1 / 3 sin 2,/3 #. (47) 

Below we show that it is true for Ti n cr 1 / 3 as well. 

First of all, excluding and from Eqn. (28) and 
using (32) and (33) we obtain the following equation 


d£ _ 2 dS 

dr tan # dr 

1 / <9t + dy~ A 1 d^ h 

2a sin # \ dr dr ) 2aX tan # dr 

It gives 

-_ ?_ 5 , A(t+ + 7~) + (7 b - Tin) cos # = 

tan # 2 Acr sin # 

+ M. 

cr sin 6 sm 6 


(48) 


(49) 


where £( 0 , #) = 1(d)/ sin# and 1(9) describes the disturbance 
of the electric current. Expression (49) corresponds to the 
conservation of the total energy flux along a magnetic field 
line. On the other hand, combining Eqns. (32)-(35), one 
can obtain the expressions corresponding to conversation of 
the ^-component of the angular momentum for all types of 
particles 


6 = 
S = 


£ ^db 7± 


,, 1 b 

£/+ 4A^ 7 



Q,r sin 6 
c 

Qr sin 0 
c 



± 4Acr 7il 


4A a 


(50) 

(51) 


Now, neglecting the difference between t + and 7 ~ so 
that £ + = £“ and T + = T _ = T (which can be done in case 
of a 3> 1 and A > 1), we obtain 


and, hence, 


<5 = ef. 


. fir sin # , 

T ( 1---) = Tin- 


Accordingly, from (51) we have 


T 


1 - 


fir sin # b 


Q = Tin- 


(52) 

(53) 

(54) 


sin # , c . 

+ ?^ ~ °- 


(55) 


On the other hand, Eqn. (33) for Act 1 gives 

_ . C 

r d6 r 

Further, using the definitions (19)-(20) to describe 7 in 
terms of £ r and and neglecting £ 0 , one can write down 

- *TTT (56) 

As a result, combining relations (53), (55), (56), and (56), 
one can finally find the following algebraic equation on the 
hydrodynamical Lorentz-factor of the secondary plasma 7 


2t 3 - 2cr 


K + ^2 + — 
2x 1 a 


T 2 + cr sin 2 # +o-= 0, (57) 
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where x = r/Rh and 

K(r,6) = 2 cos 06- sin 6^- - 7 ~ 7in cos# + ^ 77 . (58) 
d9 2a A sinfl 

Equations (57)-(58) generalize the appropriate ones 
obtained by Beskin & Rafikov (2000) for two-component 
plasma and for 7 i n -C a 1 ^ 3 . As was already stressed, these 
equations allow us to find the position of fast magnetosonic 
surface corresponding to the intersection of two roots of al¬ 
gebraic equation (57). Indeed, determining the derivative 
rdy/dr, one can obtain 


8y rdK/dr — x 2 [1 — ( 7 in/ 7 ) 2 ] 
dr ^ 0 37 — 27 i n — a(2K + x ~ 2 ) 


(59) 


As the fast magnetosonic surface is the X-point, both the 
numerator and denominator are both equal to zero there. As 
a result, evaluating rdK/dr as I\, we obtain the following 
well-known asymptotic solutions (Bogovalov 2001): 


(i) Fast rotator (y; n < a 1 ' 3 ). In this case we may neglect 
the terms containing 7i n . It gives for <5f = &{rp), etc. 


<5 F « <t~ 2 / 3 , (60) 

r F « a 1 ^ 3 sin~ 1,/3 # R L , (61) 

7f = a l ^ 3 sin 2 / 3 #, (62) 


the last expression being exact. As we see, for fast rotator 
the acceleration of the secondary plasma from 7 = 7 i n to 
7 = ct 1 / 3 takes place up to the fast magnetosonic surface. 


(ii) Slow rotator ( 7 ; n > a 1 / 3 ). In this case we obtain 


<5f - 

a 7in a \ 

(63) 

r F s= 

1/2 - 1/2 „ 
a ° 7 in -Rl, 

(64) 

7f f 

a 7in- 

(65) 


For slow rotator there is no particle acceleration within the 
fast magnetosonic surface which has a spherical shape. 


4 BEAM DAMPING AND PLASMA 
ACCELERATION 


4.1 Beam damping 


In this section we determine how the energy of the primary 
beam changes as it outflows from the pulsar magnetosphere. 
As we consider the secondary plasma in the MHD limit 
where there is no longitudinal electric field, we may neglect 
the term 88/dr. This leads to the simple determination of 
from Eqn. (33): 


£ 


r dy h 
4A a dr ’ 


( 66 ) 


where we also consider the case sin# ~ cos# ~ 1. On the 
other hand, for Act 3 > 1 one can set the quantity in brackets 
in Eqn. (30) equal to zero. After that we can rewrite Eqn. 
(31) in the following way: 


^Ke 7 b ) = 4A a 




(67) 


Now, using Eqs. (53)-(54), (56) and the definition 
1 




( C ^) 2 


2( 7 b ) 2 


( 68 ) 



Figure 1. The evolution of the beam Lorentz-factor 7 b (x) for 
different initial energies m e c 2 7 ^ (see text for more detail). 


one can finally obtain the following general equation describ¬ 
ing the beam damping: 


8_ 

dx 


idp 

dx 


= A 2 


(yP n ) 2 p + 1 7 2 x 2 p + 1 x 2 J 


+ 4 ") • ( 69 ) 


Here we use the following notations 


p = 


(7 b ) 2 - (7 b ) 2 


ini 


b \2 


A = 


4A<j 


A = 


7 -7h 


(70) 


The solutions of equation (69) for different initial energy 
m e c 2 7 n 1 of a beam (and for 7 i n = 7 = n 1 ^ 3 ) are presented on 
Figure 1. As we see, for very large initial Lorentz-factor of a 
beam 7 ^ > Act (green line), its energy remains actually con¬ 
stant. The point is that the Larmor radius tl = m e c 2 y h /eB-L 
for such a large energy of a beam is much larger than the 
light cylinder radius so the electromagnetic fields can¬ 
not disturb the radial motion of the beam. But, as we see, 
for 7 in ~ A a 1 ^ 2 , when tl « -Rl (orange line), the curve 
has oscillations corresponding to rotation in the comoving 
reference frame and beam energy slightly damps. 

On the other hand, for smaller Lorenz-factors of 
the beam 7 i n < 7 ^ < A a 1 ^ 2 (red line) the beam energy 
marginally damps to the energy of secondary plasma. Here 
again the oscillations correspond to the rotation in the co¬ 
moving reference frame. For 7 ^ = 7 = ct 1 ^ 3 (black line), 
the energy of the beam remains unchanged. Finally, if the 
initial energy of the beam is relatively close to the energy 
of the secondary plasma, the beam either damps or acceler¬ 
ates to the secondary plasma energy. Although the case of 
7 in < 7 (blue line) is impossible in real pulsars, this solution 
describes the behaviour of secondary particles: if their initial 
Lorentz-factor is less than ct^ 3 , they become accelerated to 
this energy. 

For sufficiently large initial Lorenz-factors of the beam, 
one can put p 1 and neglect the second term in (69). In 
this case Eqn. (69) becomes linear, and one could solve it 
using Green’s function: 



1 


After transition to the limit x —¥ 00 , we obtain the following 
asymptotic solution: 

p{x) ~ ^ 2 [ln(x/A) — T]. (72) 

V 7out / 
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Figure 2. Evolution of the Lorenz-factor of the beam for 
7 ^ = A a. Green solid line represents numerical solution of exact 
equation (69) and black dashed line describes analytical asymp¬ 
totic solution (72). 


Here 


7 out = lim 7 ( 2 :), (73) 

x —>-oo 

and r ~ 1-10 is a constant which weakly depends on be¬ 
havior of 'y(x) and A. For constant 7 and A > 1, this value 
is equal to the Euler-Mascheroni constant T « 0.577. 

Approximating now ln(a;/A) — T ~ 1, one can finally 
obtain the following equation describing Lorentz-factor of 
the primary particles: 

(7out ) 2 = (7in ) 2 - (4Aa/ 7 ou t ) 2 . (74) 

Hence, as was already found, large enough disturbance of 
the energy of the primary beam A 7 b ~ 7 ^ is valid only if 
the condition 7 *), < Act 1 / 2 holds. For 7 ^ > Act the energy of 
the beam remains actually unchenged. Figure 2 illustrates 
the evolution of 7 b (r) according to the numerical solution 
of Eqn. (69) in comparison with analytical solution (72). As 
we see, they are in a very good agreement. To summarize, 
the beam energy stabilizes at the following values: 


b 

Tout 

— Tout ? 

b - \ 1/2 

7in < Act ' , 

(75) 

, b 

Tout < Tout 

. b 

< Tin; 

X 1/2 ^ b ^ X 

A a < Tin < ^ cr ? 

(76) 

b 

Tout 

b 

Tin? 

7in > Act, 

(77) 

where 7 ou t 

is the final 

Lorentz-factor of the 

secondary 


plasma (see below). 

Note, that if 7 ^ > Act, the pulsar total energy losses are 
less than the beam energy. This case is thus unphysical. 


4.2 Secondary plasma acceleration 

We are now in position to evaluate the final Lorentz-factor 
of the secondary plasma 7 out . For this we look for the so¬ 
lution of equation (57) in the limit x —>• 00 . Neglecting the 
disturbance of electric current we can approximate K(r,9) 
in Eqn. (58) as 

Kao = 5oo — 7 ° U 2 • ( 78 ) 

At large distances outside the fast magnetosonic surface we 
can neglect two last terms in equation (57), which gives us 


a simple expression for 7 ou t: 

7out = 7in + ct (<5oo — ^° u * Xp™ ) ’ (79) 

According to Eqns. (60), (63), <5oo ~ inax(cr _2 ^ 3 , 7 i n CT _1 ), so 
we approximate equation (79) as 

7out = max(CT 1/3 ,7i n ) - ' 7out ^ ■ (80) 

Z A 


In what follows we consider for simplicitly the case 
7in <C ct 1,/3 , but our results could be easily generalized for 
the opposite case. 

Note, that equation (80) sets 7 ou t implicitly since 7 b ut 
depends on 7 ou t: 


_.b _,b _ (7out) 2 - (7in) 2 _ 7ta„ ^ /oi \ 

7out 7in ~ — “T-Pout (7out)- (81) 

2 7in Z 


This approximation is valid only when the beam damping is 
small, i.e., for large enough 7 in . In this limit, using equation 
(72), we finally obtain 


1/3 . 

7out = ct + 


4Act 2 
7i b n 7out ’ 


(82) 


For very large 7 ^, one can put 7out = 7in = ct 1/,s , so both the 
beam and the secondary plasma do not change their energy. 
But there is a range of the value 7 ^, where the second term 
in Eqn. (82) is larger than the first one. For such parameters 


A 1/3 2/3 

7°ut - , bu/3 ■ (83) 

(7i„) 1/3 

Finally, for sufficiently small 7 ^, one can put 7 b ut = ct 1 ^ 3 . 
Then, for A > 1 equation (80) transforms into 

7out = a 1/3 + ||. (84) 

Thus, for small enough initial energy of the beam, the energy 
of the secondary particles linearly increases with 7 ^ 1 . 

One can now construct the following approximate ex¬ 
pression, which is valid in both limits 


1/3 1 

7out = CT + 


(yiy 1 , / a i / 3 ct 2 / 3 

V2A ) + 'v(7, b J 1/3 



(85) 


and determine the maximum gamma factor of the secondary 
plasma: 


max 

Tout 


ct 1,/2 , for 7 , b = Act 


- w 1 / 2 

M • 


( 86 ) 


Figure 3 shows obtained asymptotic behavior (85). The 
exact location of the maximum is 7 b , = 6 3 ^ 4 ActJ / 2 and 
7 out x = <7 m /3 + 0.25(3/8) 3/4 ct^ /2 « ct y + 0.48ct^ / 2 . One 
should also note that the maximum gamma-factor of the 
secondary plasma is close to the initial one: even for the 
fastest pulsars, 


max / 

7out hi: 


1/6 

= CT 


10 . 


(87) 


All the expressions written above remains the same in 
case 7 in Z§> ct 1/,s as well. The only differences is in first terms 
of r.li.s of Eqns. (82), (84), and (85), which are replaced on 


1 For 7 !^ < 7 ; n , equation (80) gives 7 ou t < 7 i n , i.e., the secondary 
plasma loses energy and accelerates the primary particles 
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Figure 3. The final Lorentz-factor of the secondary plasma ac¬ 
cording to Eqn. (85). 

1 /3 

7 in instead of cfjJ . This implies that for such pulsars the 

effect of the pulsar wind acceleration via beam particles is 

1/2 

even smaller. For 7 i n 3> cr^ this effect vanishes. 


5 CONCLUSION 


It was demonstrated that the effective deceleration of the 
primary beam takes place only for relatively small beam 
energy 

7in < Acr 1/2 . (88) 


As A and a are connected with the relation (Beskin 2010) 


<j\ 



(89) 


where Wtot is the total energy losses of a pulsar and 
Wa = rrieC 5 /e 2 s=s 10 1 ' erg/s, one can easily estimate the to¬ 
tal number of pulsars having effective damping of the pri¬ 
mary beam. Setting 7 ^ = 10 7 and A = 10 4 , we obtain that 
among 2003 known radio pulsars with measured P and P, 
only 65 sources have relatively small initial Lorentz-factors 
Tin < AiT 1 / 2 . 

Thus, one can conclude that effective deceleration of the 
primary beam takes place only for sufficiently energetic pul¬ 
sars satisfying the condition Wtot > 10 36 erg/s (Crab, Vela). 
In this case during its motion through the magnetosphere 
the primary beam decelerates, accelerating simultaneously 
the secondary particles. But such an additional acceleration 
is not actually effective enough as the maximum Lorentz- 
factor of the secondary plasma does not exceed cr 1 ^ 2 , which 
is only on the factor <r 1//6 larger than the energy of the sec- 
onary plasma near the fast magnetosonic surface. For ordi¬ 
nary pulsars the deceleration of the primary beam can be 
neglected, so the energy of the primary beam as well as sec¬ 
ondary plasma remain constant. 

It is very important that the deceleration/acceleration 
takes place on the scale of the fast magnetosonic surface 
r F ~ cr 1 /, 3 f?L, he., at the distances much larger than the ra¬ 
dius of the light cylinder. This implies that even for the fast 
young radio pulsars the energetic particles of the beam will 
intersect the magnetospheric current sheet as their veloci¬ 
ties exceed the velocity of the secondary plasma which just 
determines the radial velocity of the current sheet. Indeed, 


for fast pulsars the scale tf ~ 10 2 .Rl corresponds to several 
wave length of the current sheet. 
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APPENDIX A: RADIATION FORCE DAMPING 


Inside the light cylinder the damping of a beam is defined by 
radiation losses. Outside the acceleration gap the change of 
the bulk Lorentz-factor of a beam is described by following 
equation (Zheleznyakov 1977) 


2 

me 


d7 

d7 


2 e 2 4 

3 J?| 7 ’ 


(Al) 


where R c is the curvature radius of particle trajectory (i.e., 
the curvature radius of the magnetic held line), and l is a 
path. For simplify we consider here a dipole structure of 
magnetic held inside the light cylinder. 

Dipole magnetic held lines in polar coordinates are 
described by the well-known relation r = r 0 sin 2 6, 
where ro > Rl is a parameter of the line. The appropri¬ 
ate expression in Cartesian coordinates looks like y = 
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it\/(rox 2 ) 2 / 3 — x 2 . It gives for the curvature radius 
a 1 /^4-3a 2 / 3 ) 3 / 2 


R c = r 0 - 


(6 - 3a 2 / 3 ) 


where a = x/ro- Accordingly, 


dl= r -± 


4 - 3a 2/3 


da. 


It gives 

d 7 _ 6r e f 

7 4 r 0 J 


R L /r 0 


3 V a 2 / 3 (1 - a 2 / 3 ) 

(2 - a 2/3 ) 2 


da 


R/r 0 


a (4 — 3a 2 / 3 ) 5 / 2 (l — a 2 / 3 ) 1 / 2 ' 


(A2) 


(A3) 


(A4) 


Here r e = e 2 /m e c 2 is the classical electron radius, and P is 
the stellar radius. This equation can be easily solved analyt¬ 
ically, and we obtain 


11 r e 

—--= 24 -f- 

7o Tin R L 


2 ( ^ ~f[ + 1 ) + F(r 0 , R) 


(AS) 


where F(ro,R) ~ 1 and, hence, may be neglected. As the 
initial Lorentz-factor of the primary beam after acceleration 
inside thin gap is estimated as 10 8 -10 9 (Beskin 2010), one 
can neglect the second term in the l.h.s. of Eqn. (A5), and 
we finlly obtain 


7o ~ 0.1 



10 7 . 


(A6) 


where loqj = dnriGje 2 /m e = 2f2t Ob- It gives 

dstab = ( n y /2 7 1 / 2 7 b 

Rl \UJB J A 1 / 6 

= 1500B 6 _1/2 P“ 1/2 7^7 2 “ 1/2 A6 1/6 » 1. 


(B6) 


Here we assume the value of the magnetic filed near the light 
cylinder to be of order 10 6 G which is natural for fast radio 
pulsars. Thus, one can conclude that one can neglect the 
two-fluid instability within distances which were considered 
in this paper. 

One can find even more precise expressions considering 
the values lol, w P oc r _1 as depending on the radius r, as 
it takes place in Michel solution (2)-(3). The destruction 
distance could then be found from 


Rl "b^stab 

/ 


«L 


o»p /3 (0 w b /3 (0 dr 

7(7 h ) 1 / 2 c 


= 1. 


(B7) 


Assuming now that 7 P ,b do not depend on r, we obtain 


log 


/ Rl + rfstab \ _ _ C77 b _ 

V 7?l )~ f? L ojp /3 (Rl)ooI/ 3 (Rl ) 

= 1500B 8 _1/2 P s _1/2 7772“ 1/2 ^6 V6 > !• 


(B8) 


Thus, in reality not only the condition d s tab/PL 1, but 
much more strong one 


log(dstab/ Rl'} > 1 (B9) 


is valid. 


APPENDIX B: BEAM INSTABILITY 


In this section we consider the effects leading to the beam 
instability and estimate the distance were the beam deceler¬ 
ation by the wave-particle interaction takes place. For sim¬ 
plicity, we use the following expression for plasma dielectric 
permittivity 


e(w,k) = 1 - 


K) 2 

7 3 (w - kv p ) 2 


K) 2 

( 7 b ) 3 (u; - kv b ) 2 ’ 


(Bl) 


where v p and Vb are the velocities of the secondary plasma 
and the primary beam respectively, and cu p and oib are cor¬ 
responding plasma frequencies: 


2 47re n p b 
w„,b = -■ 


(B2) 


For k || Vb || v p , the spectrum of plasma oscilla¬ 
tions (i.e., Langmuir waves) corresponds to the condition 
e(w, k) = 0. Considering now the solution as 


oj = kvb + 8u>, 


(B3) 


where 5w kvb, one can obtain for the maximum growth 
increment 


8^Jrnax — R 


K) 1 / 3 K ) 2/3 




(B4) 


with n w 1. 

Before becoming unstable, the beam travels the dis¬ 
tance 


dstab — 


8u„ 


c 7 1/2 7 b 

wGJ A 1 / 6 ’ 


(B5) 
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